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SECTION  I 


INTRODUCTION 


A.  introduction 


A  cruise  missile  is  a  vehicle  that  spends  the  major  portion  of  its 
flight  at  essentially  constant  altitude  and  speed.  Some  cruise  missiles 
are  flown  at  very  low  altitudes  in  order  to  avoid  detection.  The  ter¬ 
rain-following  capabilities  of  such  a  missile  and  its  guidance  system 
are  of  importance.  It  is  desirable  that  the  missile  should  follow 
closely  the  contour  of  the  terrain,  which  may  be  fixed  (land)  or  time- 
varying  (ocean  waves).  A  previous  study  (Ref.  1)  demonstrated  that 
optimal  control,  theory  could  usefully  be  applied  to  calculate  the  best 
achievable  accuracy  of  terrain-following.  Reference  1  showed  that,  even 
with  an  optimal  guidance  system,  terrain-following  accuracy  was  limited 
by  considerations  of  the  cruise  missile's  inertias,  airspeed,  and 
aerodynamic  configuration..  In  Ref.  1,  attention  was  focussed  on  these 
vehicle-centered  considerations,  and  little  effort  was  devoted  to  the 


and  little  effort  was  devoted  to  the 


questions  of  (i)  how  such  an  optimal  guidance  system  can  be  designed, 
and  (ii)  how  much  improvement  in  terrain-following  accuracy  can  It 
achieve  over  a  conventional  guidance  system.  The  present  repi  t 
addresses  these  questions. 


B.  CONTENTS  OF  THIS  REPORT 


Section  II  summarizes  the  results  presented  in  Ref.  1.  That  refer¬ 
ence  considered  three  alternative  aerodynamic  configurations  for 
advanced  cruise  missiles: 

1)  Conventional  wing  plus  aft  tail;  for  three  alter¬ 
native  wing  areas 

2)  Joined  wing  (a  new  type  of  wing  described  in 
Section  II);  for  three  alternative  wing  areas 

3)  Hybrid,  i.e.,  joined  wing  plus  canard;  for  one 
wing  area 
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Section  II  also  discusses  the  spectra  for  waves  and  atmospheric 
turbulence  employed  in  Ref.  1  and  in  the  present  study.  Finally,  Sec¬ 
tion  II  summarizes  the  principal  results  obtained  in  Ref.  1  with  regard 
to  the  maximum  achievable  accuracy  of  terrain-following. 

Section  III  presents  the  compensation  network  rquired  by  the  optimal 
system,  and  shows  that  it  is  stable. 

Section  IV  considers  the  effects  of  alternative  terrain  spectra, 
comparing  the  performance  of  systems  optimized  for  "rolling"  and  "rough" 
land  terrain  versus  systems  designed  to  follow  waves. 

Section  V  compares  the  performance  of  optimal  systems  versus  some 
systems  designed  by  standard  methods  (i.e.,  not  employing  optimal  con¬ 
trol  theory).  It  is  shown  that  the  latter  are  much  inferior  in  perfor¬ 
mance  to  optimal  systems. 

Section  VI  discusses  multi-controller  systems,  e.g.  systems  employ¬ 
ing  two  or  more  control  surfaces,  e.g.,  elevators  and  flaps.  It  is 
shown  that  the  airframe  characteristics  limit  the  performance  improve¬ 
ments  that  can  be  achieved  by  mu' ticontroller  systems. 

Section  VII  presents  the  conclusions  of  the  study. 
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SECTION  II 


BACKGROUND 

A.  OBJECTIVES 

The  work  reported  herein  was  performed  in  conjunction  with  a  com¬ 
panion  study  dealing  with  an  application  of  a  new  type  of  wing,  known 
as  the  joined  wing  (Ref.  1).  The  joined  wing,  as  a  general  concept, 
involves  the  combination  of  two  wings,  a  fuselage,  and  a  fin,  such  that 
the  wings  form  a  diamond  shape  both  in  plan  view  and  front  view  (see, 
for  example,  Fig.  1  overleaf). 

The  specific  application  of  Ref.  1  treats  an  advanced  cruise  missile 
capable  of  being  launche  from  a  standard  21  inch  diameter  submarine 
torpedo  tube.  The  specific  design  goals  were: 

•  Low  wave  (or  terrain)  clearance  to  avoid  detec¬ 
tion. 

•  Low  radar  signature  to  avoid  detection, 
o  Robustness  of  the  folded  vehicle. 

•  Reliable,  rapid,  and  symmetric  unfolding. 

•  Simplicity. 

•  Low  cost. 

•  High  sustained  "g"  capability  for  terminal  maneu¬ 
vers  . 

»  Long  range. 

•  Cruise  at  M  =  0.7,  little  emphasis  on  tran  onic 
and  supersonic  capability. 

This  report  deals  with  the  first  of  these  items;  it  treats  the  control 
laws  needed  to  Insure  a  low  wave  clearance  capability. 
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gure  1.  Comparison  of  Conventional  and  Joined  Wing  Configurations 


B.  BACKGROUND 


1.  The  Vehicles 

A  "cruise  missile"  is  defined  as  a  vehicle  that  spends  the  major 
portion  of  its  flight  in  a  "cruise"  mode,  i.e.,  flying  at  nearly  con¬ 
stant  altitude  and  nearly  constant  speed,  using  aerodynamic  lift  to 
support  its  weight.  The  advanced  cruise  missile  configurations  dis¬ 
cussed  in  Ref.  1  have  launch  weights,  payloads,  and  folded  dimensions 
similar  to  those  of  Tomahawk  class  vehicles.  The  specified  low-altitude 
cruise  flight  condition  (M  =  0.7  at  sea  level)  is  also  similar  to  that 
of  the  Tomahawk.  That  study  was  not  restricted  to  vehicles  comparable 
to  the  Tomahawk;  for  example,  configurations  were  considered  having 
Tomahawk-like  fuselages  but  wings  that  are  two  or  three  times  larger  in 
area  than  those  of  the  current  Tomahawk,  to  improve  maneuverability. 

In  order  to  compare  the  joined  wing  versus  the  conventional,  wing- 
plus-tail  arrangement,  3  configurations  of  each  type  were  selected.  The 
conventional  configurations  (denoted  3GA,  6GA,  9GA)  had  gross  wing  areas 
of  respectively  11. 5,  23.0,  and  34.5  ft^,  as  shown  in  Fig.  la.  The 

joined  wing  configurations  (denoted  3F,  6F,  9F)  had  gross  horizontal 
projected  areas  of  10.0,  20.0,  and  30.0  ft^.  Figure  lb  shows  the  plan- 
forms  of  these  configurations.  Configuration  3GA  is  a  baseline  configu¬ 
ration,  resembling  both,  but  not  identical  to  either  the  Vought  or 
General  Dynamics  Tomahawk  missiles.  Configurations  3GA,  6GA,  9GA  have 
aspect  ratios  ot  9.(>:  Configurations  3F,  6F ,  9F  have  aspect  ratios  of 
8.35. 

An  additional  configuration  was  also  studied.  This  configuration, 
denoted  as  (>H,  Ls  shown  in  Fig.  le.  it  is  a  hybrid  configuration  com¬ 
bining  a  joined  wing  with  a  canard. 

The  coni igurat ion  numbers  denote  very  approximately  the  nonsustained 
maneuver  capability  of  each  configuration  at  M  »  0.7,  S.L.,  maximum 

weight  conditions.  Thus,  Configurations  3GA,  IF  can  pull  approximate Ly 
3  g's,  Configurations  OF,  hGA,  bit  can  pull  approximately  b  g's,  etc. 
These  values  are  not  necessarily  sustained  g' s .  The  engine  thrust  may 
not  be  sufficient  to  overcome  the  airframe  drag  during  a  steady  maneuver 
at.  the  -it's  1  gna ted  number  of  g's. 
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The  initial  cruise  weight  of  all  configurations  is  2864  lb.  Thus 
the  more  maneuverable  configurations  sacrifice  some  fuel  because  of  the 
larger  weight  of  their  lifting  surfaces. 

2.  Airframe  Design  for  Terrain  Following 

Current  and  projected  cruise  missiles  employ  downward-looking  ter¬ 
rain-following  systems  which  measure  the  missile's  clearance  above  the 
terrain,  i.e.,  h^  in  Fig.  2.  A  bias  setting  is  incorporated  in  the 
guidance  logic,  corresponding  to  a  desired  constant  clearance  hg»  The 
difference  between  hg  and  hc,  at  any  instant,  is  the  height  error,  hg. 
If  he  has  a  symmetrical  statistical  distribution,  the  mean  value  of  hg, 
E(he),  may  be  assumed  to  be  zero.  To  achieve  stealth  hg  should  be  low, 
and  upward  deviations  from  hg  should  be  small.  With  a  low  hg,  downward 
deviations  from  hg  must  also  be  small  to  reduce  the  probability  of  hit¬ 
ting  the  terrain  (usually  called  "probability  of  clobber,"  pc).  There¬ 
fore  a  reasonable  criterion  for  the  guidance  system  is  that  it  should 
minimize  the  mean  square  height  error,  E(hg).  This  criterion  has  the 
advantage  of  being  mathematically  convenient,  and  it  is  widely  employed 

in  Optimal  Control  Theory.  The  theory  shows  how  to  compute  the  guidance 

2 

system  that  (for  a  given  airframe)  yields  minimum  E(hg). 

2 

While  there  are  advantages  to  employing  minimum  E(hg)  as  a  criterion 

2 

for  terrain  following,  it  should  be  realized  that  the  minimum  E(he)  sys¬ 
tem  is  only  an  approximation  to  the  minimum  p£  system.  This  is  because 


Figure  2.  Cruise  Missile  Flying  Above  Terrain 
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pc  depends  on  the  frequency  characteristics  of  h£  as  well  as  its  statis¬ 
tical  amplitude  distribution.  For  example,  consider  two  Gaussian  hg 
time  histories,  having  identical  E(he)  but  with  the  first  having  a 
higher  bandwidth  than  the  second.  For  equal  flight  durations,  the  high- 

bandwidth  h^  will  have  the  greater  number  of  exceedances  of  any  speci- 

o 

fled  boundary  value.  Although  E(he)  is  a  less  precise  measure  of  system 
merit  than  p  ,  it  provides  a  valid  and  cost-effective  basis  for  compari¬ 
son  of  alternative  preliminary  configuration  designs  and  was  used  for 
that  purpose  in  Ref.  1. 

3.  System  Description 

The  Ref.  1  mathematical  models  of  the  airframe  the  atmospheric 
turbulence  environment,  and  the  "terrain"  (sea  surface),  are  reviewed 
next.  The  overall  closed  loop  configuration  is  shown  in  Fig.  3. 


A1  MOSPHEHIC 
TURBULENCE 


Figure  3.  Terrain-Following  System  Block  Diagram 


a .  Airframe  Equations  of  Motion 

Short-period  longitudinal  equations  were  employed;  it  was  assumed 
that  airspeed  variations  were  suppressed  perfectly  by  an  outer  throttle 
loop.  The  equations  of  motion  were  written  in  terras  of  height  perturba¬ 
tions  (h)  and  pitch  angle  perturbations  (0)  as: 


— 

- 

- 

r—  “ 

—  - 

s-Zw 

-sh 

ZS 

<5  + 

-Zw+ZqS/U0 

-M^.s-Mw 

82-(Mq4Mi)8-M0 

9 

m6 

-MwHiqS/U0-M^s 

<5  denotes  a  single  control  surface  deflection,  and  w  is  the  vertical 

O 

component  of  turbulence.  The  terms  Z s/UQ>  M^/'U  multiplying  w  are 

*  6 
Following  customary  practice  the  Z term  is 


"gust-gradient"  terras, 
neglected,  since  «  Uq. 


The  principal  transfer  functions  of  interest  are  those  describing 
the  height  response  to  control  and  also  to  vertical  gust  components. 


-Z 6 3 2  +  Zg (Mg  +  Mj)s  +  ZgM5  -  M6Za 
-  (Mq  +  Mi  +  Zw)s  -  ^  +  MqZw] 


h_  Zw<s  -  2Mq> _ 

wg  (Mq  +  Mi  +  Zw)s  -  Mq  +  MqZw] 


(3) 


b .  Sea  Spectrum 


Reference  1  employed  a  modified  Bretschneider  wave  spectrum  appro¬ 
priate  to  a  short-crested  sea,  of  Sea  State  5,  and  a  20  knot  windspeed. 
For  a  missile  traversing  this  sea  at  M  =  0.7,  flying  into  the  inertial 
wind  direction,  the  wave  height  encountered  frequency  characteristics 
are  described  by  the  following  power  spectrum: 


whe  re 


813.9  <4 
+  6 . 3^J  ‘ 


hw  “  Wave  height  above  mean  sea  level,  ft 
u>e  *  Frequency  of  encounter,  rad/sec 


(4) 
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Figure  4  graphs  this  spectrum.  As  shown  in  Ref.  L,  it  is  a  good 
approximation  to  the  exact  Bretschneider  spectrum,  which  is  described  by 
a  more  complicated  mathematical  expression.  Note  that,  it  is  of  rela¬ 
tively  na bandwidth,  with  most  of  the  power  concentrated  in  the 
region  2  '.0.0  rad/sec.  This  is  different,  from  typical  land  cerrain 

spectra  \  i  do  not  attenuate  at  low  frequencies.  The  rms  wave  height 
is  9.13  ft 


c.  TurbuJ. mi  I'  Spectrum 


The  wo i  known  Drvden  spectrum  for  atmospheric  turbulence  was 
employed.  I  r  the  given  flight  condition,  assuming  an  rms  value  of  w^ 
of  2.9  fp  its  gives: 


■  =  (14.0438)2(u)2  +  4.5143)2 

gWg  (w2  +  7.8192)2 


(5) 


This  spectrum  is  also  shown  on  Fig.  4. 


d .  Stab ility  Derivatives 


Stability  derivatives  were  calculated  in  Ref.  1  for  each  of  the 
configurations  3GA,  6GA,  9GA,  6F,  9F,  and  6H  shown  in  Fig.  1.  In  all 
cases  a  positive  static  margin  was  assuci  i  as  listed  in  Table  1,  which 
also  gives  weights  and  inertias.  Table  2  lists  the  derivatives  and  the 
h/5  and  h/w  transfer  functions  in  the  forms: 

o 


h 


h 


Ah5s  +  Bh6s  +  ch6 
s2^2  +  Bs  +  C) 


A.  s  +  B. 

_  s _ s 

s2  (s2  +  Bs  +  C) 


(6) 


For  more  background  on  the  Bretschneider  wave  spectrum  see  Ref.  2. 
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TABLE  1 


FLIGHT  CONDITIONS  FOR  TERRAIN-FOLLOWING  ANALYSIS 


Mach  No. 
Altitude 


0.7 

For  all  configurations 

Sea  I.evel 


-  50.5% ,  for  Configurations  3GA,  6GA,  OCA,  3F,  6F,  9F 

%  Total  Fuel  Used 

-  6.2X,  for  Configuration  611 


Conf  igura  i  ion 

3GA 

3F 

6GA 

6F 

6H 

9GA 

9F 

Weight  (lb) 

2234.9 

2077.69 

2396.8 

2183.52 

2771  .5 

2622.7 

2320.58 

Pitch  Inertia 

1620.0 

1591 .03 

1815.6 

1644.2 

2363.54 

1983.5 

1789.7 

Static  Margin  (ins.) 

4.175 

2.012 

5.907 

1.992 

2.0 

7.23 

2.0562 

Reference  Area*  (ft^) 

11.5 

10.0 

23.0 

20.0 

15.9992 

34.5 

30.0 

CL  at  1  g 

0.2673 

0.2858 

0.1434 

0.1502 

0.238 

0.1046 

0.1064 

g  Per  Degree  of 

Control  Deflect  ion' 

1.1341 

1.1907 

1.5478 

3.0099 

2.5817 

2.8766 

4.5531 

*The  reference  area  is  the  area  employed  to  define  C^.  For  Configurations  3GA,  6CA,  9GA 
it  equals  the  gross  wing  area,  for  Configurations  3F,  6F,  6H,  and  9F  it  equals  the  gross 
horizontal  protected  area  of  the  joined  wings  (excluding  the  canard  for  611). 

^"Control"  «  all-moving  tail  for  30A,  6CA,  9GA;  wing-warping  (measured  at  front  wing  root) 
for  3F,  6F,  9F,  and  all-moving  canard  only  for  611. 


TABLE  2,  DERIVATIVES  AND  TRANSFER  FUNCTIONS 

(Flight  Conditions  Specified  in  Table  1) 


CONFIG.  3GA 

CONFIG.  6CA 

CONFIG,  9CA 

CONFIG.  3F 

CONFIG.  6F 

CONFIG.  9F 

CONFIG.  611* 

MALPHA 

-14.30900000 

-36.14000000 

-60.76000000 

-3.97450000 

-7.36460000 

-10.00340000 

-5.24880000 

MW 

-0.01830000 

-0.04620000 

-0.07770000 

-0.00510000 

-0.00940000 

-0.01280000 

-0.00671200 

ZALPHA 

-959.8800000 

-1,790.900000 

-2,345.450000 

-538.7485000 

-1,003.024900 

-1  ,444.244800 

-864 .7700000 

ZW 

-1.22760000 

-2.29040000 

-3.14040000 

-0.75300000 

-1  .38510000 

-1 .84710000 

-1 . 10600000 

MA  DOT 

-0.12480000 

-0.22230000 

-0.29620000 

-0.07560000 

-0.21080000 

-0.41330000 

-0.27900000 

MW  DOT 

-0.00015860 

-0.00028430 

-0.00037880 

-0.00009670 

-0.00026960 

-0.00052850 

-0.00035700 

MDKLTA 

0. 59050000 

1 .12900000 

1.56460000 

0.25000000 

0.70510000 

1  .21540000 

0.57070000 

Z DELTA 

1.81860000 

3.81859000 

5.60510000 

-3.13910000 

-6.03610000 

-8.43150000 

-1  .31800000 

MQ 

-0.40800060 

-0.72450000 

-0.98000000 

-0.25180000 

-0.70270000 

-1.37750000 

-0.69800000 

D 

1.76040000 

3.23720000 

4.41660000 

1 .00040000 

7.29860000 

3.63790000 

2.08300000 

C 

14 .80986000 

37.79939480 

63.83759200 

4 .16410540 

8.33790977 

12.54778025 

6.01998800 

A,  H  DELTA 

-1 .81860000 

-3.81859000 

--5. 60510000 

3.13910000 

6.03610000 

8.43150000 

1 .31800000 

B,  H  DELTA 

-0.96895000 

-3.61544101 

-7.1  532  2862 

1.02774134 

5.51397735 

1 5.09913020 

1  .28768600 

C,  H  DELTA 

540.7867926 

1,863.922257 

3,501.231194 

159.6634780 

808,0943191 

1 ,839.678797 

500.4411030 

A,  II  UC 

-1 . 22760000 

-2.29040000 

-3 .1404  0  0  00 

-0.75300000 

-1 .38510000 

-1  .84710000 

-1.1 0600000 

B,  H  WC 

-1 .001721  60 

-3.31878960 

-6.1  5518400 

-0.37921080 

-1  .94  661954 

-5.08876050 

-1  .  54397600 

*Wlth  canard  1/4-chord  sweep  angle  -  35  deg. 
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Figure  4.  Command  and  Disturbance  Spectra 

C.  LIMITS  ON  OPTIMUM  TERRAIN-FOLLOWING 
PERFORMANCE  (SINGLE  CONTROLLER) 

For  a  given  airframe,  the  guidance  system  can  be  designed  to  mini- 
2 

mize  E(hg)  through  the  use  of  optimal  control  theory.  This  yields  the 

2 

guidance  system  that  gives  the  minimum  achievable  E(h£)  for  a  specified 

2 

level  of  mean  square  control  deflection  E(6  ). 

2 

It  might  be  thought  that  any  desired  level  of  E(hg),  even  zero, 

2 

could  be  obtained  by  allowing  E(6  )  to  increase.  This  is  not  true  for 
conventional  (aft  tail)  configurations  in  which  elevator  is  employed  to 
control  height.  Such  configurations  have  height-to-elevator  transfer 
functions  which  incorporate  one  or  more  right-half-plane  zeros.  (Trans¬ 
fer  functions  of  this  kind  are  called  "nonmiuimum  phase.")  For  such  a 
system  there  is  a  finite  minimum  achievable  mean  square  error  (Ref.  1). 


ft 
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2.  Optimum  Single-Controller  System 


Reference  1  derives  the  equations  for  the  optimal  single-controller 

"M 

system,  i.e.,  the  system  that  minimizes  the  performance  index,  J. 

J  =  E(h2)  +  k2E(62)  (7) 

2  2 
For  a  giver,  k  we  wish  to  know  the  value  of  E(he)  obtained,  and  the 

2 

associated  mean  square  control  deflection  E ( 6  ).  As  shown  in  Ref.  1, 
these  values  are  given  by  expressions  of  the  type: 

E(.heJopt  =  2tTJ  f .  +  w2<i>ddw2}  ds  (g) 


i  r ■5°°  -  - 

Associated  E(62)  =  j  {wa$RRWa  +  Wb<|)ddWb}  ds 


In  each  expression  the  first  term  gives  the  contribution  to  J  that  stems 
from  the  portion  of  he  or  6  that  is  correlated  with  the  terrain.  The 
second  term  gives  the  contribution  due  to  the  portion  of  he  or  6  that  is 
correlated  with  the  atmospheric  turbulence.  Only  two  terms  appear 
in  each  expression:  there  are  no  cross-correlation  terras  involving 

^Rd’  ^dR  because  it  has  been  assumed  that  the  terrain  (waves)  and 
turbulence  are  uncorrelated. 

E(h?)  and  the  associated  E(<$2)  were  calculated  in  Ref.  1  for  each 
c  opt 

of  the  configurations  in  Fig.  1.  The  results  are  presented  below. 


3.  Conventional  (Aft  Tail)  Configurations 

2 

Figure  5  shows  the  mean  square  height  error,  E(he),  for  aft  tail 
configurations  3GA,  6GA,  and  9GA  graphed  versus  the  corresponding  mean 

2 

square  control  deflection  E(6  )  required  by  the  optimal  control  system. 


12 
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Figure  5.  Optimal  Terrain-Following  Performance  of 
Single-Controller  Aft  Tail  Configurations 


Values  of  the  weighting  parameter  k  are  indicated  on  the  graphs.  (The 

control  deflection.  6,  indicates  deflection  of  the  all-moving  tail 

surfaces  in  units  of  degrees.)  At  high  values  of  k  (l„e.,  "expensive 

control")  corresponding  to  the  left  side  of  the  graphs  the  control 

deflections  are  small  and  have  little  effect  on  the  missile,  which  flies 

more  or  less  in  a  straight  line  with  the  preponderant  part  of  the  height 

2 

error  being  due  to  the  waves.  Thus  the  asymptote  of  E(he)  for  k  +  <*>  '  s 

2 

essentially  equal  to  the  mean  square  wave  height,  81.2  ft  .  As  s 

decreased  the  system  accuracy  increases,  but  for  reasonable  valut  >f 
2  o 

E ( 6  ),  e.g.,  of  the  order  of  10  deg4,  the  gain  in  a.  curacy  is  slight. 

2  2 

For  example,  with  F,(6  )  =  10  deg  Configuration  3GA  reduces  its  mean 
square  error  only  to  74/81.2  =  0.91  of  its  open-loop  value.  If  much 
larger  E(6  )  could  be  produced  cr  if  the  control  effectiveness  param¬ 
eters  could  be  increased  the  accuracy  could  be  improved,  but 

even  with  infinite  mean  square  control  deflections  the  best  accuracy 
that  can  be  achieved  with  airframe  3GA  is  only  E(hg)  =  22  ft^,  i.e., 
0.27  of  the  open-loop  mean  square  error. 
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This  situation  is  only  slightly  improved  by  selecting  a  larger  wing 
area,  as  shown  by  the  graphs  for  Configurations  6GA  and  9GA  in  Fig.  7. 
Nor  is  any  significant  improvement  gained  by  reducing  static  margin. 
This  expedient  was  studied  by  recalculating  the  optimum  system  per¬ 
formance  with  Ma  set  equal  to  zero,  and  the  remaining  derivatives 

,  o 

unchanged.  For  any  given  k  the  zero  static  margin  values  of  E(hg)  and 

9 

E(<5  )  were,  when  graphed  on  Fig.  5,  virtually  coincident  with  the  points 
graphed  for  the  normal  static  margins  listed  in  Table  1. 

It  is  clear  from  the  above  results  that,  even  with  ideal  guidance 
system  dynamics  (e.g.,  no  actuator  lags,  no  digital  processing  delays), 
the  mean  square  control  deflections  required  for  close  terrain-following 
of  the  assumed  wave  spectrum  are  prohibitively  large  for  the  aft  tail 
configurations . 

4.  Atmospheric  Turbulence  Response 

The  Dryden  spectrum  employed  to  describe  low  altitude  turbulence 
(Eq.  5)  is  widely  employed  in  aircraft  guidance  and  control  studies. 
The  mean  square  height  error  due  to  turbulence  was  much  smaller  than 
that  due  to  terrain  for  all  the  configurations  studied.  Therefore  the 
effect  of  varying  the  Dryden  spectrum  is  secondary,  hence  no  parametric 
variations  of  atmospheric  turbulence  characteristics  were  performed  In 
Ref.  1,  or  in  the  present  study. 

5.  Terrain-Following  Performance  of  Single- 
Controller  Joined  Wing  Configuration 

The  joined  wing  configurations  3F,  6F,  9F  obtain  pitch  control  via 
wing-warping,  a3  described  in  Section  3  of  Ref.  1.  With  a  parabolic 
variation  of  twist,  corresponding  to  uniform  torsional  modulus  GJ  along 
both  wings,  the  center  of  pressure  of  the  lift  induced  by  twist  Is  ahead 
of  the  c.g.,  and  hence  the  height/twist  transfer  function  Is  minimum 
phase.  However,  the  pitching  moment  arm  is  small  for  Configuration  3F. 
For  the  hybrid  configuration,  6H  (Fig.  3c),  a  long  moment  arm  Is 
obtained  by  employing  an  all-moving  canard  as  the  primary  pitch  control. 
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Figures  6  and  7  shows  the  terrain-f ollowirg  performance  of  these 

configurations  when  optimally  controlled.  The.  gust  and  wave  spectra  are 

identical  to  those  employed  previously.  Comparing  Figs.  6  and  7  with 

Fig,  5  it  is  seen  that  the  joined  wing  and  hybrid  configurations  yield  a 

mean  square  height  error  which  tends  to  zero  as  k  +  0.  Unfortunately, 

as  was  the  case  with  the  aft  tail  configuration,  the  magnitudes  of  the 

2 

control  deflection  E(6  )  are  impracticably  large  for  close  terrain¬ 
following  of  the  assumed  short-crested  sea.  The  reason  for  this  is 
explained  below. 

Consider  Eq.  10,  which  gives  the  h/6  transfer  function  for  Configu¬ 
ration  6H,  where  6  =  deflection  of  the  all-moving  canard,  measured  in 
degrees : 


h  =  1.318s  +  1.287687s  +  500.441103 

5  s2(s2  +  2.083s  +  6.019988) 


h_  _  1 .318(s  +  0.4885  ±  19.4797  j) 

6  s  2  (  s  2  +  2.083s  +  6.019988) 


(10) 


(ID 


The  Bede  diagram  of  Eq.  10,  given  on  Fig.  8,  exhibits  a  pronounced  dip 
around  19.5  rad/sec,  due  to  the  light  damping  of  the  numerator  fac¬ 
tors.  This  dip  attenuates  the  h/6  response  in  the  region  where  it  is 
most  needed,  i.e.,  in  the  region  around  the  peak  of  the  wave  spectrum. 
If  the  dip  is  eliminated  (or  smoothed)  the  response  improves  drama¬ 
tically.  This  can  be  demonstrated  by  calculating  the  mean  square  height 
error  and  mean  square  control  deflections  for  a  smoothed  transfer  func¬ 
tion  : 


yielding  E(hg) 


h 

6 


100 

O 


45.1  ft2  and  E(_62j  =  6.8  deg2. 


(12) 
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Figure  8.  Hetght/Canard  Deflection  Transfer  Functions 
for  Hybrid  Configurations 

The  results  presented  up  to  this  point  of  the  present  report  were 
previously  given  in  Ref.  1.  That  reference  is  primarily  concerned  with 
airframe-centered  aspects  of  Cruise  Missile  Design,  and  does  not  inves¬ 
tigate  the  practical  realizability  of  the  Theoretical  Optimum  Guidance 
Systems.  By  contrast,  the  present  report  is  not  solely  concerned  with 
the  optimum  level  of  performance,  but  also  addresses  the  question  of 
whether  the  optimum  compensation  can  be  constructed  from  practical  com¬ 
ponents.  Tiie  next  section  shows  that,  although  the  theoretical  optimum 
compensation  cannot  be  achieved  in  practice,  it  can  be  closely  approxi¬ 
mated,  with  negligible  performance  penalty  being  Induced  by  the  approxi- 
rant ion . 
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Phase 


SECTION  III 
COMPENSATION  NETWORK 


A.  INTRODUCTION 

The  compensation  network  required  by  the  optimal  controller  is  shown 
to  require  ideal  lead.  The  governing  equations,  described  in  Appendix  A 
of  Ref.  1,  are  first  reviewed. 

B.  SIMPLIFIED  SYSTEM  MODEL 

A  simplified  model,  which  ignores  the  gust  disturbance  input,  is 
given  in  Fig.  9.  R  represents  the  "sea  input",  Gc  is  the  compensation 
network  and  G  represents  the  open  loop  dynamics  of  h/ 6 .  From  Fig.  9, 
write 

6  =  [1  +  GCG]_1  GCR  =  WaR  (lb) 


Figure  9.  Simplified  System  Model 


and  recall  that  represents  the  quantity  to  be  determined  in  the  opti¬ 
mization  process.  Specifically,  the  Wlener-Hopf  equation  of  interest  is 


1  +  —  GG 
k2 


Wa^RR  ,  •;  ^RR  “ 
k~ 


(17) 


TR- 1147-2 


18 


For  the  modified  Bretschneider  spectrum,  the  solution  is  obtained  by 
letting  G  =  N/D  where  N  is  2nd  order,  D  is  fourth  order) 


Wa 


D 

A 


__N _ s2 _ 

A^2  (s+6,3)4 


(s+6.3)4 


(18) 


D  represents  the  open  loop  roots  and  A  represents  the  optimal  closed 
loops.  From  Appendix  A,  Ref.  1, 


Ns- 


A(s+6,3)4  (s+6.3)4 


where  L  is  a  cubic  polynomial.  Therefore; 

cl 


DLa 

W  =  _ 

a  k2As2 


(19) 


(20) 


For  the  assumed  cruise  missile,  IV  is  7th  order  over  6th  order. 

a 

Next,  solve 


wa  =  h  +  GcGr'G 


(21) 


for  Gc : 


DL 


Wa  I  1 


GW, 


k2  (As2  -  NLa) 


(22) 


We  see  that  the  numerator  ot  C  in  7  th  order  while  the  denominator  is 
6th  order,  indicating  pure  lead  at  high  frequency.  The  practical  sig¬ 
nificance  of  this  requirement  can  be  assessed  by  considering  the  bode 
diagrams  of  G£  for  two  typical  examples,  as  discussed  below. 
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C.  SPECIFIC  EXAMPLES 


The  Bode  plots  for  Gc,  are  shown  in  Figs.  10  and  11  for  the  aft-tail 
Configuration  9F  and  the  joined  wing  plus  canard  Configuration  6H, 
respectively.  Each  Gc  is  stable,  as  determined  by  a  Routh-Kurwitz  Test, 
i.e.  all  the  poles  lie  in  the  left-half-plane. 

Within  the  frequency  region  1  <  m  <  50  rad/sec  the  Bode  diagrams  can 
be  closely  approximated  by  Bode  diagrams  of  systems  without  pure  lead, 
i.e-,  systems  having  transfer  function  numerators  of  equal  or  lower 
order  than  the  corresponding  denominators.  Outside  this  frequency 
region  the  amplitude  of  the  wave  spectrum  is  small,  hence  the  practical 
inability  to  generate  the  pure  lead  demanded  by  the  optimum  Gc  is  not 
important.  Note  also  that  for  w  >  50  rad/sec  an  accurate  solution  re¬ 
quires  a  more  complex  airframe  and  control  system  model  than  considered 
here.  Such  a  model  would  include  actuator  lags,  digital  processing  lags 
and  aeroelastic  modes.  If  such  high-frequency  effects  were  included  the 
optimum  Gc  would  change.  Nevertheless  it  would  always  be  possible  to 
avoid  Gc  having  a  higher  order  numerator  than  denominator  by  appropriate 
adjustment  to  the  high  frequency  characteristics  of  the  wave  spectrum. 
Such  adjustment  would  involve  a  negligible  loss  of  performance. 

From  the  above  discussion  we  reach  the  following  conclusions  regard¬ 
ing  the  optimum  compensation: 

1)  It  Is  stable. 

2)  It  cannot  be  physically  realized  exactly,  due  to 
the  requirement  for  a  pure  lead  (i.e.,  an  ideal 
differentiator). 

3)  it  can  be  closely  approximated  by  a  physically 
realizable  compensation  with  negligible  loss  of 
performance. 

D.  SECTION  SUMMARY 

The  optimal  compensation  Is  a  ratio  of  polynomials  in  the  complex 
frequency  variables.  The  numerator  is  seventh  order  whereas  the  denom¬ 
inator  is  sixth  order,  indicating  the  need  for  approximation  of  the 
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Compensation  for  9F  Cor.fi; 


optimal  compensation  at  high  frequencies.  Bode  plots  for  the  compensa¬ 
tion  networks  were  given  for  the  9F  and  6H  configurations.  Each  network 


SECTION  IV 


EFFECT  OF  ALTERNATIVE  TERRAIN 


A.  INTRODUCTION 

A  short-crested  sea  spectrum  was  selected  as  the  "terrain"  to  be 
followed  in  the  examples  presented  in  the  previous  sections.  It  is 
clear  from  Figs.  5,  6  and  7  that,  even  for  optimal  systems,  this  terrain 
is  hard  to  follow.  Large  control  deflections  are  required  to  obtain 
close  terrain-following.  For  the  nonminimum  phase  (conventional  wing- 
plus-tail)  configurations,  even  infinitely  large  control  deflections 
cannot  reduce  the  mean  square  height  error,  E(h  ) ,  much  below  20  ft  , 
i.e.,  approximately  25  percent  of  the  open-loop  mean  square  error.  In 
this  section  we  examine  the  effect  of  alternative  terrain  spectra. 
Reference  1  briefly  discussed  the  asymptotic  behavior  of  the  nonminimum 
phase  airframes  following  various  types  of  terrain,  i.e.,  their  behavior 
with  unlimited  control  deflections.  Reference  1  did  not  study  how  the 
accuracy  of  terrain-following  obtainable  with  finite  control  deflec¬ 
tions  is  affected  by  the  terrain  spectrum.  This  important  question  is 
addressed  below. 


B.  ALTERNATIVE  TERRAIN  SPECTRA 


Various  land  terrain  spectra  are  in  use  (see  Ref.  3).  The  simplest 
of  these  is  the  first-order  form: 


Urr)+ 


/la [ E ( R)2] 
s  +  a 


(23) 


o 

In  this  equation  E(R  )  is  the  mean  Bquare  terrain  altitude  variation  and 
a  is  the  encounter  break  frequency,  which  depends  on  the  terrain  rough¬ 
ness  and  missile  speed  as  follows: 


a 


(24) 
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where 


Vj  =  Missile  Inertial  Speed  (FPS) 

ag  =  Spatial  Correlation  Distance  (FT),  i.e.,  the 
average  distance  between  two  sample  terrain 
points  which  differ  in  magnitude  by  a  factor 
of  1/e  =  0.368. 

Here  we  shall  consider  two  types  of  terrain: 

1)  Rolling  terrain  with  a  =  2,000  ft 

s 

2)  Rough  terrain  with  ag  =  500  ft 

For  both  terrains  the  mean  square  height  variation  was  chosen  as 
2 

81.2  ft  .  This  value  equals  the  mean  square  wave  height  variation  of 
the  wave  spectrum  analysed  in  previous  sections.  It  is  a  convenient 
choice  because  it  removes  the  effect  of  terrain  amplitude  when  making 
comparisons  between  the  results  to  be  presented  in  this  section  and 
those  given  in  earlier  sections.  Since  the  optimum  system  is  linear, 
the  results  may,  of  course,  be  scaled  for  any  terrain  amplitude  of 
interest.  The  terrain  spectra  then  become: 

For  Rough  Terrain 

i»«kJ+  '  rrraia  <25> 


For  Rolling  Terrain 


UrrJ+ 


7.9584 
s  +  0.39 


(26) 


C.  TERRAIN-FOLLOWING  PERFORMANCE 

Figure  12  presents  the.  results  of  Wiener-Hopf  optimizations  for  Con¬ 
figurations  6GA  and  GF,  following  rolling  terrain.  The  abscissa  is  the 

2 

expected  value  of  the  mean  square  control  deflection  E(<5  ).  Care  must 

2 

be  used  in  interpreting  E(6  )  since  the  control  surfaces  employed  for 
the  joined  wing  Configuration  6F  are  different  from  those  used  in  the 
conventional  Configuration  6GA.  For  Configuration  6F  6  denotes  the 


number  of  degrees  of  wing  warp  applied  at  the  wing  root.  A  reasonable 
limiting  value  based  on  maximum  "g"  structural  limits  for  Configura¬ 
tion  6F  is  -15  <  6  <  15  degrees.  For  Configuration  6GA  6  denotes  the 
deflection  of  the  horizontal  tail,  measured  in  degrees.  Typically,  for 
similar  maximum  "g"  this  6  would  be  in  the  range  -5  <  6  <  5  degrees. 

7  2 

Thus,  reasonable  maximum  levels  of  E(5“)  would  be  100  deg  for  Con- 

2 

figuration  6F  and  and  11  deg  for  Configuration  6GA.  The  correspond- 

2  2 

ing  levels  of  terrain-following  accuracy  are  E(he)  «  17  ft4  for  6GA  and 
2 

E(heJ  -•  11.8  ft^  for  Configuration  5F. 

Comparing  Figs.  5,  6,  and  12  shows  that  the  rolling  terrain  can  be 

followed  much  more  closely  than  the  short-crested  sea,  for  a  given  mean 

2 

square  control  deflection.  In  particular,  for  low  values  of  E(6  )  of 
the  order  of  1  deg,  the  rolling  terrain  permits  a  substantial  improve¬ 
ment  over  open-loop  constant-altitude  flight,  whereas  for  the  short- 
crested  sea  only  a  slight  improvement  is  obtained.  As  shown  in  Fig.  13, 
for  rough  terrain  the  results  are  intermediate  between  the  results 
obtained  for  waves  and  for  rolling  terrain. 

The  optimum  terrain-following  performance  of  the  larger  crnise  mis¬ 
siles  (Configurations  9GA  and  9F)  is  shown  in  Figs.  14  and  15.  The 

2 

smaller  configurations  (3GA  and  3F)  yield  optimum  E[heJ  and  E[6^J  as 
illustrated  in  Figs.  16  and  17.  Different  definitions  of  6  for  the 
joined  wing,  hybrid  and  conventional  configurations  makes  it  difficult 
to  perform  meaningful  comparisons  between  the  configurations.  Neverthe¬ 
less,  it  is  clear  that,  if  adequate  control  power  is  available,  the 
minimum  phase  Configurations  3F,  6F,  9F,  6H  have  superior  performance. 
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Figure  14.  Rolling 


Terrain;  9G A.,  9F 


Figure  13.  Rough  Terrain;  9CA,  9F 
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SECTION  V 


COMPARISON  OF  OPTIMAL  AND  CONVENTIONAL  CONTROL 
SYSTEM  PERFORMANCE 

A.  INTRODUCTION 

Section  II  of  this  report  has  shown  that  none  of  the  optimally 

controlled  configurations  considered  achieved  close  following  of  the 

short-crested  sea  without  large  control  deflections.  By  definition,  the 

2 

optimal  control  system  yields  a  smaller  E[heJ  than  any  other  system.  It 
is  therefore  of  interest  to  compare  the  optimal  system  performance  as 
given  in  Figs.  5,  6  and  7  with  the  performance  of  conventionally 

designed  control  systems,  i.e.,  systems  which  are  designed  by  standard 
systems  synthesis  procedures  (as  opposed  to  optimal  control  theory). 
Our  objective  is  to  determine  how  much  is  gained  by  employing  an  optimal 
instead  of  a  conventional  guidance  system. 

For  simplicity,  instead  of  the  exact  airframe  transfer  functions 
given  in  Section  H  we  shall  assume  that  the  h/6  transfer  function  is 
of  the  form  100/s".  This  is  a  reasonable  approximation  to  the 
exact  transfer  functions  except  that  it  omits  the  Bode  amplitude  "dip" 
around  20  rad/sec.  As  shown  in  Section  II,  elimination  of  this  dip 
improves  the  performance  of  the  optimally  controlled  system,  yield- 
lng  E[heJ  =6.8  deg  with  h/ 5  =  100/s  .  These  numbers  should  be  borne 
in  mind  when  reading  the  results  presented  below. 

Two  "Conventional"  designs  will  be  studied:  first,  a  washout  compen¬ 
sator,  second  an  ideal  lead.  The  former  can  be  physically  realized 
exactly,  the  latter  only  approximately,  so  the  results  given  for  the 
ideal  lead  compensation  are  slightly  optimistic,  but  nevertheless 
approximately  indicate  the  level  of  performance  that  could  be  attained 
by  a  conventional  system. 
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B.  WASHOUT  COMPENSATOR 


Consider  a  washout  network  in  the  feedforward  path,  as  defined  by 

Fig.  18.  In  Fig,.  18,  R  represents  the  terrain  input,  E  is  the  attitude 

error'  and  h  is  the  altitude.  Table  3  tabulates  E(6^),  E(h^)  and  A(s) 

2  2 

for  the  closed  loop  system,  and  Fig.  19  graphs  E(h  )  versus  E(<$  ). 

Figure  19  shows  that  the  use  of  a  washout  network  leads  to  a  peak  in 
the  plot  of  E(h^)  vs.  E(<5^). 

Table  4  compares  the  washout  design  against  an  optimally  controlled 

2 

system  with  the  same  transfer  function  (100/s  ). 

The  results  given  on  Table  3  (and  Table  4)  assume  a  A(s)  representa¬ 
tive  of  optimally  designed  systems,  i.e.,  critically  damped.  For  a  pole 

2  2 

assignment  of  A  =  (s+3)  +  (5)  ,  an  underdamped  set,  the  overshoot  is 
higher;  E(S)^  =  9.75  deg^  ,  E(h^)  =  130  ft^.  It  is  clear  that  the  wash¬ 
out  compensated  conventional  design  is  considerably  inferior  to  the 
optimally  designed  system. 

C.  IDEAL  LEAD  COMPENSATION 

It  was  shown  in  Section  III  that  the  optimal  design  requires  lead 

compensation.  Therefore,  it  is  worth  checking  effects  of  an  ideal  dif¬ 

ferentiator  as  a  compensator  for  a  conventional  design  (refer  to 
Fig.  20).  The  values  of  E(h^)  and  E(5^),  as  well  as  A(s),  the  closed 

loop  poles,  are  tabulated  in  Table  5  and  graphed  in  Fig.  21.  It  can 

2  2 

be  seen  that  the  plot  of  E(h  )  vs.  E(<$  )  monotonical ly  decreases  as 

2 

E(6  )  increases.  This  is  apparently  a  property  of  lead  compensation. 

2 

By  way  of  comparison,  the  optimal  system  requires  2,528  deg  of  mean 


Figure  18.  W,,  ;hout  Compensation  Feedback  System 
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TABLE  3.  CONVENTIONAL  DESIGN  USING  WASHOUT 


E(S2) 

E(h2) 

A 

0.0000033 

83.89 

(s+0.2)2  +  (O.l)2 

0.03328 

91.135 

(s+1)2  +  (l)2 

0.505 

104. 

(s+2)2  +  (2)2 

2.288 

112.1 

(s+3)2  +  (3)2 

6.622 

112.85 

( s+4 ) 2  +  (4)2 

12.787 

108. 

(s+5)2  +  (5)2 

83.524 

69.315 

(s+10)2  +  CIO)2 

331.0523 

28.3 

(s+20)2  +  (20)2 

628.915 

14.6 

(s+30)2  +  ( 30 ) 2 

1593.543 

4.17 

(s+60) 2  +  (60) 2 

3750.38 

1 .01 

(s+125)2+(125)2 

TABLE  4.  COMPARISON  OF  OPTIMAL  DESIGN  WITH  WASHOUT 


DESIGN 

E(«2) 

E(h2) 

A 

Optimal 

0.0369 

74.9 

(s+2.24)2  +  ( 2 . 24  )2 

Washout 

0.7774 

106.73 

(s+2.24)2  +  (2 .24  )2 

Optimal 

1.475 

54.5 

(s+5)2  +  (5)2 

Washout 

12.787 

108.0 

(s+5)2  +  ( 5 ) 2 
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Figure  20.  Ideal  Lead  Compensation 


TABLE  5.  IDEAL  LEAD  COMPENSATION 


E(S2) 

E(h2) 

A 

0.0166 

83.8 

s  +  0.01 

6.529 

77.16 

s  +  2 

37.873 

59.304 

s  +  5 

129.01 

37.3 

(3  +  10) 

29  39.  12 

1.37 

(s  +  100) 

1791.73 

0.9 

(.3  +  125.7) 
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Figure  21.  System  Performance  with  Modified  Bretscbneider , 
Rolling  Terrain,  and  Pure  Lead  Compensation 


square  control  deflection  to  produce  0.9  ft"'  of  E(h2).  The  "pure  lead" 
"conventional  design  requires  3,791.73  deg  -  Again,  the  evidence  is 
clear;  the  short-crested  sea  described  by  the  modified  Brentschneider 
spectrum  is  hard  to  follow  at  M  -  0.7,  even  with  ideal  lead  compensa¬ 
tion,  and  such  compensation  requires  larger  control  deflections  than  are 
required  by  the  optimal  system. 

A  physical  explanation  of  the  peak,  in  the  plot  of  E(h  )  vs.  E(<3  ) 
tor  the  conventional  design  is  readily  apparent.  For  E(6  )  =  0  the  con¬ 
troller  is  in  a  mode  where  it  follows  an  inertial  altitude  reference. 

2 

For  small  values  of  E(6  ),  with  no  lead,  the  vehicle  responds  out  of 

2 

phase  with  the  input  (loop  delay  is  too  large)  and  E(h  )  actually  rises 
above  that  for  following  an  inertial  reference.  As  E ( 6  )  increases,  the 
bandwidth  of  the  closed-loop  system  becomes  sufficiently  wide  so  that 
ECh1")  drops  below  the  Inertial  reference  value. 

In  this  regard  (see  Table  3)  note  that:  the  bandwidth  of  the  conven- 

2 

tional  design,  for  a  given  E^hgj,  Is  much  larger  than  the  optimal 

2  2 

design.  For  example,  for  a  E^hgj  on  the  order  of  70  ft.  the  bandwidth 
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of  the  washout  design  is  of  the  order  of  10  rad/sec,  whereas  the  optimal 
design's  bandwidtn  is  only  2  rad/sec. 

I 

y  y 

The  tendency  for  the  plot  of  E(h£j  vs.  E(6  )  to  peak  is  greatly 
diminished  if,  instead  of  the  short-crested  sea  spectrum,  the  spectrum 

y  y 

for  rolling  terrain  is  employed.  Figure  19  compares  E(h|J,  E(f>  )  for 
the  two  spectra.  This  confirms  that  the  peaking  occurs  because  of  the 
added  phase  lag  incurred  when  no  lead  compensation  is  used. 

D.  SUMMARY 

|i  2 

Conventional  compensation  produces  a  peak  in  the  plot  of  E(h  )  vs. 

y 

E(<5  ).  Lead  compensation  produces  a  monotonically  decreasing  relation¬ 
ship.  The  performance  of  the  two  conventional  designs  is  notably 
inferior  to  the  system  designed  by  optimal  control  theory. 

¥ 


* 


» 
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SECTION  VI 


SYSTEM  SURVEYS  FOR  MULT ICONTROLLE R  SYSTEMS 

A.  INTRODUCTION 

The  term  "Multicoatroller"  will  be  used  here  to  denote  systems 
employing  two  or  more  separate  controls,  e.g.,  flaps  plus  elevator.  In 
Ref.  1  it  was  suggested  that  good  terrain-following  performance  could  be 
obtained  by  controlling  an  "inner  loop"  by  means  of  flaps  such  that  the 
effective  "outer  loop"  transfer  function  controlled  by  the  elevator  had 
a  desirable  form  (e.g.,  K/s  )  for  which  the  h/6  Bode  diagram  did  not 
have  a  notch  near  the  wave  spectrum  peak,  frequency.  Reference  1  did  not 
explore  this  possibility  in  detail,  although  some  discussion  was  given 
of  how  a  height/elevator  transfer  function  of  the  form  100/s  could  be 
achieved  by  employing  flaps  on  Configurations  6H  and  6GA.  It  was 
concluded  that  the  hybrid  Configuration  6H  could  readily  attain  the 
desired  inner-loop  transfer  function,  whereas  this  was  more  difficult 
for  the  conventional  6GA  cor  figuration. 

Although  a  multicontroller  system  can  offer  improved  performance 
over  a  single-controller  system,  the  former  is  necessarily  more  complex, 
perhaps  prohibitively  so  for  cruise  missiles  which  must  satisfy  severe 
constraints  on  cost  and  complexity.  The  decision  as  to  whether  or  not 
to  employ  a  multicontroller  system  therefore  depends  on  whether  it  can 
provide  a  considerable  Increase  in  terrain-following  accuracy.  This 
question  is  studied  below.  The  performance  of  optimally  controlled 
systems  with  h/6  transfer  functions  of  various  forms  is  assessed,  and 
desirable  forms  are  indicated.  The  results  of  this  section  are  of 
value  in  deciding  whether  che  additional  cost  and  complexity  of  a  multi¬ 
controller  system  is  justified  by  the  increased  performance  that  it 
provides . 

B.  SURVEY  OF  K/s 

The  mean  square  height  errors  and  mean  square  control  deflections 
for  K/s  open-lo^p  dynamics  and  the  short-crested  sea  spectrum  are 
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tabulated  in  Table  6.  Also  listed  is  the  closed-loop  optimal  pole. 
From  Table  6  is  is  seen  that  an  open-loop  sys  .em  with  K/s  dynamics  can 
produce  a  small  mean  square  height  error  (about  10  deg  of  control 
deflection  with  about  3  ft2  error).  Moreover,  the  closed-loop  bandwidth 
of  30  rad/sec  constitutes  a  reasonable  set  of  dynamics  for  an  "outer" 
altitude  loop. 

C.  SURVEY  OF  K/s2 

9 

The  survey  of  K/s  (Table  7)  indicates  that  this  form  jf  transfer 
function  is  less  desirable  than  K/s.  The  control  effecciven  terms 
must  be  extremely  high  (K  of  the  order  of  4000)  to  achieve  an  error  on 
the  order  of  1  ft  .  Moreover,  the  closed-loop  natural  frequency  of 
approximately  65  rad/sec  is  somewhat  high  for  an  altitude  outer  loop. 


TABLE  6.  SURVEY  OF  K/s 


K 

E(_62J 

Elh2) 

A 

10 

15.96 

20.1 

s  +  10 

20 

13.78 

8.59 

s  +  20 

30 

9.6 

3.25 

s  f  30 

40 

6.71 

1.46 

s  +  40 

100 

1.5 

0.074 

s  +  100 

TA3LE  7.  SURVEY  OF  K/s2 


K 

E  (  6  2  J 

E[h2] 

A 

10 

3.69 

74.9 

(s  +  2. 24) 2  +  ( 2 . 24 ) 2 

20 

4.912 

65.85 

(s  +  3.2)2  +  (3 . 2  )2 

50 

5.9 

54.5 

(s  +  5)2  +  (5)2 

100 

6.8 

45.1 

(s  +  7. 1)2  +  ( 7 . 1 ) 2 

1000 

6.2 

6.13 

(s  +  22)2  +  (22)2 

4000 

1.58 

0.9 

(s  +  45)2  +  (45) 2 
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D.  SURVEY  OF  K/s3  AND  K/s4 


For  K/s3,  a  value  of  K  =  103  is  required  to  produce  E(_62)  =  3.3 
and  E[h“)  =  5.8  ft  .  The  closed  loop  roots  are 

A  =  (s  +  46.4) [(s  +  23. 2)2  +  (40. 2)2]  (27) 

For  K/s4,  a  value  of  K  =  10®  is  required  to  produce  E(62J  =  0.8, 
E[h2 J  =  1.8  ft2  and 

A  =  [(s  +  38. 3)2  +  (92.4)2][(s  +  92. 4)2  +  (38. 3)3]  (28) 


E.  IMPLICATIONS  OF  SYSTEM  SURVEYS  FOR 

AIRFRAME  SELECTION 

The  K/s  transfer  function  is  indicated  to  be  the  most  desirable 
form,  since  good  performance  is  attained  for  a  reasonably  low  value  of 
K.  We  now  examine  the  practical  realizability  of  such  a  transfer 
function.  As  shown  in  Fig.  22  it  is  assumed  that,  in  addition  to  the 

outer  loop  control  <$i,  an  additional  control  $2  is  employed..  In 

Fig.  22  62  is  denoted  fs  "Flap",  but  62  may  in  fact  be  a  combination  of 

several  aerodynamic  surface  controls,  e.g.,  for  Configuration  6H  a 

combination  of  front  wing  flaps  and  rear  wing  flaps.  From  Table  5  we 
put  K=3Q  as  a  desirable  value.  Then,  following  the  procedure  of  Ref.  1, 
p.  112,  the  crossfeed  transfer  function  Cx  is  given  by  the  solution  of 


(Mq  +  M;js  - 


[Ma  -  (,^62^252l +  ^hg  ^ 
s^g2  +  Bs  +  Cj 


where  B,  C,  ,  B^^  ,  etc-  are  given  in  Table  1. 
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For  Conf iguration  6H  this  yields 


C 


30s3  +  30Bs2  +  30 C  -  A^s2  -  B^s  - 

Z62{S2  -  (,Mq  +  MqJ s  -  [ Mq  -  iM62/Z62JzcJ} 


(30) 


30s3  +  61.172s2  +  179.3119s  -  500.441103 
z62{s2  +  °*977s  +  5.2488  -  864.77  (_M62/Z§2J} 


(31) 


For  the  crossfeed  to  be  constructed  physically  C  must  not  be  unstable. 
This  leads  to  the  same  condition  on  as  derived  in  Ref.  1,  i.e., 

Ms2  „  5.248 

Z  5  2  <  *5^ 


or,  approximately 


^62 

z62 


0 


which  implies  that  the  flap  deflection  which  generates  an  increase  in 
lift  must  also  generate  a  nose-up  pitching  moment,  or  at  least  zero 
pitching  moment.  In  other  words,  a  "Direct  Lift"  type  of  control  is 
necessary.  This  is  more  readily  achieved  with  a  joined  wing  or  canard 
configuration  than  with  a  conventional  aft  tail  configuration. 

A  second  condition  that  must  be  met  by  the  crossfeed  transfer 
function  is  that  C  should  be  physically  realizable.  This  implies  that 
the  numerator  of  C  should  be  of  equal  or  lower  order  than  the 
denominator.  Clearly,  C,  as  calculated  above  does  not  meet  this 
condition.  It  is  therefore  necessary  to  modify  C  by  introducing 
additional  high-frequency  denominator  terms.  These  additional  terras 
will  not  seriously  degrade  the  system  performance,  because  the  wave 
input  peak  frequency  is  relatively  low  (6,3  rad/sec). 
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SECTION  VII 


CONCLUSIONS 


1.  Optimal  control  theory  of  the  Wiener-Hopf  type  can  conveniently  be 
applied  to  the  design  of  cruise  missile  terrain-following  guidance 
systems.  For  single-controller  systems  the  optimal  compensation 
transfer  function  is  stable,  but  may  have  a  numerator  of  higher 
order  than  its  denominator.  Thus,  it  must  be  approximated  if  it  is 
to  be  physically  constructed.  The  performance  loss  due  to  this 
approximation  is  small. 

2.  For  multi-controller  systems  to  achieve  good  performance  at  least 
one  of  the  airframe  control  surfaces  should  provide  either  zero  or 
nose-up  pitching  moment  when  the  control  is  deflected  so  as  to 
increase  lift.  This  is  more  easily  achieved  by  a  joined  wing  or 
canard  configuration  than  by  a  conventional  aft  tail  configuration. 

3.  Of  the  three  types  of  terrain  considered  (rolling  terrain,  rough 
terrain,  and  a  short-crested  sea)  the  short-crested  sea  provided  the 
most  severe  operating  environment,  and  required  large  control 
deflections  for  close  terrain-following.  The  effect  of  atmospheric 
turbulence  was  relatively  small. 

4.  Terrain-following  systems  designed  by  optimal  control  theory  give 
appreciably  better  performance  than  those  obtained  by  standard 
system  design  methods.  When  following  the  short-crested  sea  the 
latter  systems  give,  mean  square  errors  greater  than  would  have  been 
obtained  by  flying  at  constant  inertial  altitude.  The  control  sys¬ 
tem  feedback  laws  for  existing  cruise  missiles  were  not  available  to 
the  authors,  so  it  is  net  known  how  existing  cruise  missile  guidance 
systems  compare  with  the  standard  (non-optimal)  systems  analysed  in 
this  report.  If  existing  cruise  missile  guidance  systems  employ  any 
of  the  standard  (non-optimal)  feedback  laws  analysed  in  this  report 
considerable  improvements  in  terrain-following  accuracy  would  be 
obtained  by  changing  to  optimal  feedback  laws. 
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